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Abstract 

In this paper, we introduce further our recently designed family of L 2 orthogonal Space-Time codes for CPM. 
<^ With their advantage of maintaining both the constant envelope properties of CPM, the diversity of Space-Time 

codes and moreover orthogonality, and thus reduced decoding complexity, these codes are also full rate, even for 
more than two transmitting antennas. 

The issue of power efficiency for these codes is first dealt with by proving that the inherent increase in 
bandwidth in these systems is quite moderate. It is then detailed how the initial state of the code influences the 
coding gain and has to be optimized. For the two and three antennas case, we determine the optimal values by 
computer simulations and show how the coding gain and therewith the bit error performance are significantly 
improved by this optimization. Q 

I. Introduction 

t/5 ! The combination of space-time coding (STC) and continuous phase modulation (CPM) is an attractive 
field of research because the combination of STC and CPM brings the benefits of diversity and low 
power consumption to wireless communications. Zhang and Fitz [1] were the first to apply this idea by 
^. . constructing a trellis based scheme. This idea was pursued by Zajic and Stiiber in [|2|| for full response CPM, 
00 \ afterward optimized and extended to partial response CPM in But for these codes the decoding 
^ | effort grows exponentially with the number of transmitting antennas. This was circumvented by burst- 
■ wise orthogonality introduced by Silvester, Schober and Lampe in J5]| and by block-wise orthogonality 
£Ni '• established by Wang and Xia [6JL7J. Unfortunately, based on the Alamouti code [8], this latter design is 
restricted to two antennas. This restriction was partially circumvented in |[9l by using quasi orthogonal 
STC for 4 transmitting antennas. However, to our knowledge, there was so far no consistent orthogonal 
design for a general number of antennas. 

In our previous publications, by relaxing the orthogonality condition, we have been able to construct 
^ \ new L? -orthogonal space-time codes which achieve full rate and full diversity with low decoding effort 
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c3 ' More precisely, in IfTOfl we have generalized the two-antenna code proposed by Wang and Xia from 
pointwise to L 2 -orthogonality and in [11 J and Ifl2ll . we introduce the first L 2 -orthogonal code for CPM 
with three antennas. In this paper we shortly display some of these results and focus on several properties 
of these codes. Of special interest is the optimization of the bit error rate which depends on the initial 
phase of the system. Our simulation results illustrate the systemic behavior of these conditions. 

II. Parallel Coded CPM 

A. Space-Time code block 

In MIMO systems, the signal sent s(t) is transmitted via L t transmitting antennas. Here, similarly to 
0, the continuous phase modulated sending signal will be split into blocks 



«l,l(t) ■•• s l,L t (t) 
Si t ,l(t) ■•• s L . L Jt) 



(1) 
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where each element of the sending matrix 



S m ,r(t) = y 6XP C? 27r ^m,r(*)) ( 2 ) 

is defined within the time interval (L t l+r— 1)T <t< (L t l+r)T. The indexes m and r denote respectively 
the sending antenna and the time slot in a block. E s is the symbol energy per transmit antenna Tx and 
per symbol time T. The phase signal 

7 

<Pm,r{t) = O m {Ltl + r) + hY J d { ^}q{t - i'T) + c m , r (t) (3) 

i=i 

is essentially the phase of a conventional CPM signal [13] with an additional correction factor c mjf .(t). 
The data symbols dm$ are uniquely determined for each element of the code block. A central point of 
the coding scheme is the mapping between these symbols and the data sequence dj with dj E = 
{— M + 1, — M + 3, . . . , M — 1}. As in conventional CPM, h = 2m /p is the modulation index where 
m and p are relative primes. The phase smoothing function q(t) is a continuous function with q{t) = 
for t < and q(t) = 1/2 for t > 7T and 7 is the memory length. The phase memory will be obtained 
from the definitions of correction factor and the mapping of the data symbols. 



B. L 2 — orthogonality 

To achieve simplified blockwise decoding as in systems based on orthogonal space-time block codes 
(OSTBC) HQ, it is sufficient OH to force each block to be L 2 -orthogonal 

(l+l)L t T 

S{t)S"{t)dt = E S I (4) 



lL t T 



where I is the identity matrix. Hence the correlation between two different Tx antennas s mr (b) an d s m r r (b) 
is zeroed over the / th complete STC block if 

(l+l)L t T 

S m ,r(*)C'.r(*) d* = (5) 
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with m m'. Now, from theses conditions, L 2 -orthogonal codes can be constructed IfTOl ffTTTl ffT2l by 
expressing our design criteria: 

1) the mapping of the data symbols dm]l to the data sequence dj and 

2) the correction factor c m r (t). 

Mapping: The most convenient mapping [fTTI is the parallel mapping where each data symbol of one 
time slot r for m = 1, . . . , is mapped to the same symbol of the data sequence dj. The following L t 
data symbols at r + 1 are mapped to dj + i and so on. Each row of the code is therewith modulated by the 
same symbols of the data sequence and 

^m,r — dL t l+r-i+l (6) 

Due to this parallel structure, this solution is named parallel codes (PC). This mapping is illustrated in 
Figure [Q 

Correction factor. The correction factor is determined by the phase difference A0 between the different 
transmitting antennas at the end of each symbol in the code block. To ensure phase continuity, this 
difference has to be 

. 1/2 at t = (21 + 1)T for 2 Tx antennas 
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Fig. 1. Mapping of the uncoded data sequence to the code block for 3 Tx antennas 



. and 1/3 or 2/3 at t = (31 + 1)T and 2/3 or 4/3 at t = (3/ + 2)T for 3 Tx antennas. 
We define the linear correction factor (linPC) as 

m-1 t-(L t l + r)T 

CmA 1 ) = — t f <7) 

and a more complex one based on the phase smoothing function q(t) as 

m — 1 7 

CmAt) = ^^^r- J2 Mt - i'T). (8) 

Lt i=l 

Without loss of generality, the correction factor of the first antenna is set to Ci r (t) = 0. Using Eq. ([8]) with 
Eq. ([3]) one can see that the sums from each equation can be merged. So, by combining the correction 
factor with the data symbol, we get new alphabets for each transmitting antenna m: 

Q dm = {-A/+1+M _m+3 + M ...,M_ 1+ M}. 

Consequently, this code is called offset PC (offPC) and may be seen as L t conventional CPM signals 
with different alphabet sets VL dm for each antenna m. The alphabet of the first antenna Vt dl is equal to 
the alphabet of a conventional CPM Vt d . For example, for two transmitting antennas we have Ci r (t) = 0, 
c 2 , r = Y2=i V(t ~ ^T)/ 2 and tnen =V d = {—M + 1, -M + 3, . . . , M - 1} and Q d2 = {—M + 1 + 
1/h, —M + 3 + 1/ h, . . . ,M — 1 + 1/h}. This intuitive representation greatly simplifies modulation and 
demodulation [fTTI . 

With these definitions of the coding scheme, we can rewrite the correction factor in the more generic 
form 

9 m (L t l + r + 1) = 9 m (L t l + r)+ £(L t l + r) (9) 

where the function £(L t l+r) guarantees the continuity of the phase for any correction factor and mapping. 
For parallel mapping (similarly to conventional CPM) and c m>r (t) = 0, we get £(L t l +r) — \d 
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III. Initial Phase 

Our model includes now all the necessary parameters to construct L 2 -orthogonal STC. The modulation 
index h and the phase smoothing function q(t) can be chosen with the usual restrictions of conventional 
CPM detailed above. 

As proved by our simulations in section [TV] the values of the initial phases 9 m (l), which are known to 
have no influence in conventional CPM systems [fT3l . will be shown to have instead a great importance 
on the performance of the proposed code. 



-5 -4 -3 -2 -1 1 2 3 4 5 
normalized frequency fT 



Fig. 2. Power spectral density of the linPC analyzed with the Welch algorithm 

A. Continuous-time Model 

First, we need to introduce another formalism than the block structure used for design. The signals sent 
by each Tx antenna m are rewritten as 



s m (t) = Jj± exp (j27r[0 m (l) + h ^ d t q(t - (i - 1)T) + c m {t)}). (10) 
' * i=i 

In Eq. (flOl) the phase memory terms 9 m (L t l + r) get included in the summation term. Only the initial 
phase 9 m (l) remains. This is due to the property of continuity in the definition of every symbol over the 
whole time N C T where iV c is the number of transmitted symbols. As a result of the parallel mapping, the 
data symbols di are also equal for each antenna. For offPC codes, we use the modified alphabets detailed 
in section III-BI and no additional correction factor is necessary. For linPC codes, the correction factor 
c m (t) simplifies to a continuous linear function c m (t) = (m— l)t/L t . 

It is interesting to see that for the second Tx antenna, the correction factor causes a constant phase 
offset of 2%/L t per symbol and of 2n per block. For the 3 Tx antennas case, these offsets are multiplied 
by 2. The same effect is observed from the offset added to the alphabet VLd m . 

This phase offset induces a frequency shift. For a phase shift of 27r(m — l)/L t on a period of T, we 
get a frequency shift of Af m — (m — l)/TL t for antenna m and a symbol length T. 

Figure |2] shows the simulated power spectral density for the linPC code with 3 Tx antennas, h = 1/2, 
g = 2 and M = 4. To achieve an attenuation in power of -30dB, a bandwidth expansion of some 5 HzT is 
necessary. This corresponds to an increased bandwidth demand of ^rijf- = ^ = 0.133. As the alphabet 
size M grows, the absolute bandwidth of the CPM signal widens but the shift caused by the coding 
scheme is constant. Therewith the additional relative bandwidth eventually decreases. 

It should be recalled here that the highest achievable rate for linear codes with 4 Tx antennas is 3/4 
lfT4~l . It means that to transmit the same quantity of data during the same duration of time as for our 
proposed code a 25% increase in bandwidth is necessary. 

With this continuous-time formalism, the signals sent can be rewritten in vector form as s(t) = 
[£i(£)>£2(*)) ' ' ' ;lL t (^)] T - Below underlined variables denote the vector representation and non-underlined 
variables the previously used matrix form. Thus we can write s(t) as the product of two matrices and 
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C_(t) and a vector d(t) 



s(t) 



Es 
L t T 



QC(t)d(t) 



The matrix of initial values 6 = diag(#) and the matrix of correction factors C_(t) 
L t x L t diagonal matrices obtained from the vectors 



•exp(j27r^(l)) 
exp(j27T0 2 (l)) 

exp(j27r0 Lt (l))J 



c(t) 



exp(j'27rci(t)) 
exp(j27rc 2 (l)) 

.exp(j27rc Lt (l)). 



(11) 

diag(c(t)) are 
(12) 



As a result of the parallel mapping, the vector of data symbols can be written as 



d(t)=exp(j27rh^2d iq (t-(i-l)T))[l 1 ... l] H . 



(13) 



B. Code Performance 

The performance of theses codes may now by evaluated using the classical pair-wise error probability 
(PWEP). We assume optimal demodulation, i.e. maximum likelihood (ML) sequence detection (MLSD). 
Furthermore, it is considered that for < t < N C T the signal s(t) modulated by the data sequence dj 
is the one truly sent. The PWEP is then the probability that this signal is erroneously detected as signal 
S(t) modulated by dj 

|A T A(t)|| 



PWEP = P(s(t) -> 1(0 1 A) = Q 



(14) 



where A is the channel matrix which is assumed to to have frequency flat quasi-static Rayleigh fading and 
mutual independent elements. The energy of noise is given by N and Q(-) is the cumulative distribution 
function of the normal distribution (Q-function). The normalized difference vector A(i) is given by 



A(t) 



rAi(t)] 
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(15) 



The PWEP is minimized by maximizing the product of the eigenvalues of the signal matrix 031, 



N C T NcT N C T 

f |Ai(i)| 2 dt /Ai(i)A3(t)dt ... /Ai(t)At (t)dt 



N C T N C T N C T 

fA 2 (t)Al(t)dt f |A 2 (t)| 2 dt - /Aa(t)AJ t (*)«tt 



N C T N C T N C T 

fA Lt (t)Al(t)dt fA Lt (t)A*(t)dt ••• / |A it (t)| 2 dt 





(16) 



Using Eq. (fTTI) and (fT5l ) the signal matrix can be written as 



(a) offset PC with REC (b) offset PC with RC (c) linear PC with REC 

Fig. 3. Simulated BER for varying initial phase #i(l) and and 83(1) = for 3 Tx antennas at 13dB SNR 



N C T 



C s = / A(t)A H (t)dt 



N C T 



[QC(t)(d(t) - d(t))] [QC{t){d{t) -d{t))] H dt 



e( / c(t) Ad(t) Ad"(t)c"(t) dt)e" 



(17) 



with Ad(t) = d{t) — d(t) and the integral over the matrix acting element-wise. Since d(t) has equal 
elements and is multiplied by its Hermitian transpose, we get an L t x L t all-ones matrix. This matrix is 
multiplied with the matrices of the correction factor and we obtain the correlation matrix of the correction 
vector 



c s = ej c(t)c H (t)dt e H . 



(18) 



By writing c m (t) = H ^ i c(t) in Eq. © and ©, we get 



L t Lt 

The autocorrelation c m {t)(f m (t) is therewith always one and we have 

cxp(-^i27r£(t)) exp(--±j2nc(t)) ... exp(-^j27rc(i)) 
exp(-^-i27rc(t)) exp(-^j27rc(i)) ... exp (-- ^j27rc(f)) 



c(t)c H (t) 



cxp 



(iizlj 27 ra(t)) exp(^i27rc(t)) ... exp(-iLi2^c(t)) 



(19) 



(20) 



The elementwise integration of this matrix is easily computed [12J in the special cases where c(t) is a 
linear function (linPC codes) or a sum of raised cosines (OffPC codes) [fT2ll . In both cases, the matrix C s 
is shown to have full rank [fTTTl and thus our codes achieve full diversity. However, the PWEP approach 
doesn't provide here any valuable estimate of the coding gain. For that reason, we detail hereafter some 
statistical estimates to show the influence of the initial phase values upon the coding gain. 
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Fig. 4. simulated BER for varying initial phase 62 for 2 Tx antennas at SNR=12.5dB 
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TABLE I 

Initial phase values for minimal BER (SNR=12.5dB; 2REC) 




IV. Simulations 

In this section, we benchmark by simulations the performance of the proposed codes. In all the 
simulations, we used an alphabet size of M = 4 with 2-bit Gray-coding, a modulation index of h = 1/2, 
a memory length of 7 = 2 and 12 samples per symbol. The signal is disturbed by complex code-block- 
wise Rayleigh fading of variance one. In this section all given phase values are relative values, e.g. 1 
corresponds to 27r or 360°. 



A. Two transmit antennas 

The bit error rate (BER - E^/Nq) for the proposed two transmitting antenna codes depends on the 
difference of the initial phase A9 = 02(1) — 0i(l)- Figure @] shows the results of computer simulations 
for linPC and offPC codes with different phase smoothing functions q(t). The variation of performance 
covers almost one decade. This shows the importance of a carefully chosen initial phase. 

Mainly, the position of the minimal BER seems to depend on the correction factor used. Between linPC 
and offPC, the minima are shifted by 1/4. However, the phase smoothing function used has only minor 
influence on the position of the minima. It is also interesting to see that the distance between the minima 
is 1/2 and further simulations show a periodicity of 1. 



B. Three transmit antennas 

For the three antennas codes, the BER also depends on the initial phase m (l). Fig. [3] shows the 
simulation results for different codes with varying initial phases for the first and second antenna and null 
initial phase for the third antenna. It can be seen that the phase offset for a minimal BER depends on the 
correction factor chosen (Fig. |3(a)| and |3(c)[ ). However, similarly to the two antenna code, the form of the 
phase smoothing function q(t) has almost no influence on the position of the minima (Fig. |3(a)] and |3(bj] ). 
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Fig. 5. BER for Rayleigh-fading channel 

In Table U the optimal initial phase are summarized for the first and second antenna. Other simulations 
for different alphabet sizes M, modulation indexes h and memory length 7 validate the position of the 
minima and prove it is an important issue for the optimal design of parallel codes. 

C. Bit Error Rate 

Fig. [5] shows the influence of Rayleigh fading channels on a CPM transmitter with a different number 
of antennas. 

For the optimal two antenna system we used in the optimal case a frequency offset of A9 = 0.19 for 
linPC and of A9 = 0.4 for offPC. For the optimal three-antennas-codes, we took the following values 
from Fig. [3] 

. offPC: 0i(l) = 0.1, 2 (1) = 0.45; 

. linPC: 0!(1) = 0.4, 2 (1) = 0.15. 
The non-optimal codes have no phase offset (6>i(l) = 6 2 (1) = A6> = 0). 

The optimized codes achieve the expected performance gain. For high SNR the BER decreases with 
5dB/dec similar to a two antenna system with full diversity. The three Tx antennas code achieves a decay 
of some 3.5dB/dec. This validates the property of full diversity. 

Fig. \5\ shows clearly the improvement of the coding gain by using an optimized initial phase. Comparing 
the optimized codes with the non-optimal ones, we achieve an additional coding gain of around 5dB for 
the two antenna system and of around 7dB for the three Tx antennas. 

V. Conclusion 

In this paper, we detail the construction and analyze some of the properties of L2-orthogonal STC-CPM 
for two and three transmitting antennas. These codes are attractive due to their low-effort-decoding and 
the few restrictions the code-family set upon the parameters of CPM. We give a general formulation for 
two and three antenna parallel codes and introduce a continuous-time representation of the CPM signals. 
With this representation we are able to analyze how the coding gain depends on the initial phase of the 
system. Furthermore, we give the optimal values for the initial states obtained from computer simulation. 
The significant gain in performance for typical Rayleigh fading channels is shown and compared with 
non-optimal parallel codes. 
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